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ABSTRACT 


Reduced velocity distributions are derived from three-dimensional 
measurements of the velocity distribution of electrons in the 7-500 eV range 
in the electron foreshock. Bump-on-tail reduced distributions are presented 
for the first time at the foreshock boundary consistent with Filbert and 
Kellogg's proposed time-of- flight mechanism for generating the electron beams. 
In a significant number of boundary crossings, bump-on-tail reduced 
distributions were found in consecutive 3 sec measurements made 9 sec apart. 

We conclude that, although the beams are linearly unstable to plasma waves 
according to the Penrose criterion, they persist on a time scale of 3-15 sec. 
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Introduction 

Upstream from the earth f s bow shock, enhanced fluxes of energetic 
electrons are observed when the local magnetic field connects to the shock 
surface [Ogilvie et al«, 1971; Feldman et al . , 1973; Anderson et al., 1979; 
Feldman et al . f 1982]. This is the electron component of the earths 
foreshock [Greenstadt, 1976] and occupies the space behind the surface of 
tangent magnetic field lines as shown in Figure 1. Closely associated with the 
enhanced electron fluxes is electrostatic plasma noise that is most intense at 
the upstream edge of the foreshock [Scarf et al . , 1971; Filbert and Kellogg, 
1979; Anderson et al., 1981]. Filbert and Kellogg [1979] proposed that the 
plasma noise could be generated by the bump-on-tail class of reduced electron 
velocity distributions, although at the time, such velocity distributions had 
not yet been observed. The basis for this proposal was that the upstream edge 
of the foreshock is spatially dispersed with respect to parallel velocity due 
to the solar wind electric field drift of the particles. Electrons 
originating on connecting field lines at the bow shock with the highest 
parallel velocity move upstream along a resultant path that is closest to the 
tangent field line, while electrons with lower parallel velocity cross-field 
drift for a longer time and are carried further downstream (see Figure 1). 

This time-o f-fl ight separation produces a spatial variation with respect to 
parallel velocity and was suggested by Filbert and Kellogg as the mechanism 
for the existence of peaked reduced velocity distr lbutions at the foreshock 
boundary. Recently, Fung et al . [1 982] have suggested that the foreshock 
plasma noise can be produced by electron velocity i iscr ibutions which are not 
of the bump-on-tail class postulated by Filbert and Kellogg. The motivation 
for the distribution suggested by Fung et al . was the lack of observational 
evidence for the bump-on-tail . In this letter, we present new experimental 
findings which bear on these theoretical discussions. 

Recent election plasma measurements confirm the velocity dispersive 
boundary model; i.e., when the magnetic field configuration changes in such a 
way that the foreshock boundary is swept past the observing spacecraft, the 
energy at which the onset of upstream electrons occurs is found to decrease 
with increasing penetration into the foreshock [Anderson et al . , 1979; Feldman 
et al., 1983]. Peaks m fast (3 second) two-dimensional measurements at the 
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foreshock boundary have been reported by Feldman et al . [1983], and peaks in 
energy spectra of electrons moving upstream from the bowshock have been 
measured by Gurnett and Frank [1975]* These measurements , however, do not 
address the presence of multi-peaked reduced distributions as postulated by 
Filbert and Kellogg. It is the purpose of this paper to present the required 
fast, three-dimensional measurements and to report the first observations of 
bump-on-tail reduced distributions at the foreshock boundary. A later 
publication will present more observations in a morphological study of the 
electron foreshock. 


Observational Result 

The data presented here were obtained from the Goddard Vector Electron 
Spectrometer experiment on ISEE-1 [Ogilvie et al . , 1978]* The instrunent is 
fully three dimensional and consists of six detectors, each with a narrow 
field of view (8.5 degrees x 11 degrees), which sample the energy range 7-500 
eV in 16 approximately logarithmically spaced energy steps in 0.5 second. The 
energy cycle is repeated six times per spacecraft rotation and constitutes the 
basic measurement which is repeated each 9 seconds. In one three second 
rotation period, 576 vector measurements of the velocity distribution 
f( v x ,v y ,v z 9 1) are macle * The data are transformed to the cylindrical coor- 
dinate system which is symmetric about the measured magnetic field vector, B, 
and which is at rest in the frame moving with the electron bulk velocity. 

This reduces the four-space ,v z ) to three dimensions (f,w x ,w lt ), where 

w^ and w„ are the electron proper frame velocities perpendicular and parallel 
to B, respectively. In this way it is possible to compute the reduced 
distribution , F(w„), from measurements of the three-dimensional distribution 
using 


F(w„) = 2ir/f(w A ,w n ) w t dw f . (1) 

An example of an electron velocity distribution measured just inside the 
leading edge of the foreshock is shown in Figure 2 (upper panel) and the 
reduced distribution computed from it using Equation 1 is shown m Figure 2 
(center panel). The integration of f(w 1 ,w tt ) was performed in strips at 

7 

constant w„ spaced aw„ = 5 x 10 cm/s apart. The integration limits of were 
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extended by a smooth extrapolation to w l f : 0 as 4 co . The magnetic field 
points to the right m Figure 2, which as shown in the sketch of the geometry 
for this event in Figure 1, was directed toward the point of tangency from the 
spacecraft. Therefore, electrons moving away from the shock along the 
connected field line will have negative values of w !t on the plot. We see in 

this example that there is a peak in the contour map and a bump on the tail of f 

8 8 
F(w h ) at w lf = -7 x 10 cm/s. There are also closed contours at w lf = + 7 x 10 

cm/s which outline a local minimum xn the f(w lt w lf ) surface. There is a 

plateau in F(w n ) corresponding to this minimum which appears to mirror the 

bump on tail. This presence of a return flux of electrons toward the shock is 

not an uncommon feature at the foreshock boundary. 

Before interpreting the contours further, we briefly review the analysis 
method by which they are produced. The underlying points are the center 
coordinates of 12 degree pitch angle-speed bins into which measurements made 
during one spin have been sorted, compatible with the field of view of the 
detectors. The surface of f(w JL ,w n ) is represented in three dimensions by a 
mesh of interconnecting triangles, formed by clusters of three nearest 
neighboring data points. The f(w x ,w fl ) surface in the space between measured 
data points is linearly interpolated using the three data points at the 
triangle vertices that bound the given region. The exception to this approach 
is in the cone-like data void about the w !t -axis. This region is filled with 
data by an interpolating spline that is computed over pitch angle on constant 
speed shells, followed by triangulation of the interpolated points m the ^ame 
way as the other data. Because of the spline interpolation, the representa- 
tion of the f(w t ,w Jt ) surface along the w Tf -axis will depend on 1) the pitch 
angle gradients near the edge of the data void, and on 2) the assumed 
cylindrical symmetry about the w ff -axi s. If there are weak pitch angle gradi- 
ents, then the contours will tend to be circular as is the case for speeds 
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below 3 x 10 cm/s. On the other hand, the peak m the contour map on the 

negative w„-axis is the result of a well-developed trend in the measured data 
and the assumed symmetry. The smooth variation in the data near the edge of 
the data void extends over a large enough phase space element to include 
measurements in several 0.5 second energy sampling cycles, which suggests that 
these data are not seriously tune-aliased. 



Although the peak in this f(w A ,w„) contour map lies on the negative 
w„-axis in a region sampled only by interpolated data, the relevant velocity 
space producing the bump on the F( w„) distribution lies off the axis and 
includes a range of w A sampled by actual measurements. Vie show this in the 
bottom panel of Figure 2 in which iso-contours of the integrand of Equation 1, 
w A f(w A ,w„»), normalized for each w„» to the maximum, ( w A f(w <t ,w„») ) max , are 
plotted on the w A -w rt grid. Because of the w A weighting factor, the integrand 
has a maximum off the w A s 0 axis (contour level 1.0) and decreases with 
respect to w A on either side of the maximum (contour levels 0.8, 0.6, 0.4, and 
0.2). The contour separation broadens and extends well into the velocity 
space sampled by actual measurements in the parallel velocity strip producing 
the positive slope of the bump on F(w„) (highlighted by the two vertical 
dashed lines in Figure 2). Note also that there is a sharply defined parallel 
velocity at which this broadening occurs. The nearly vertical alignment of 
the contours means that there is a lower limit to the parallel velocity of the 
upstream electrons, w , which is a confirmation of the time of flight 
mechanism. Compare this feature with the more gradual spreading of the 
contours m the right half plane m which the solar wind heat flux electrons 
have no sharply defined lower limit in w tl . 

A second example of a bump on the reduced distribution is shown m Figure 
3. This is a case in which there was excellent data coverage close to the 
magnetic field direction. The elongated, elevated contour at the parallel 
velocity, w fe , in the f(w 1P w„) map and the w A f contour map in the bottom panel 
clearly show that the perturbations to the distribution function have a wide 
range in w A and a lower bound, w c , in w„ . This organization of upstream 
velocities according to w„ is a basic prediction of the time of flight model 
sketchei in Figure 1. 

The vertical bars on the F(w tI ) plots in Figure 2 and Figure 3 are the 685 
confidence intervals of the reduced distribution. The propagation of errors 
in determining F(w„) is complicated by the interdependence of neighboring 
elements of the f(w A ,w„) surface and oy the fact that F(w |( ) is an integral 
over that surface (Equation 1). Therefore, a Monte Carlo technique has been 
used wmch takes into account the propagation of both systematic and 
statistical errors. The expected uncertainties in electron counts based on 
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Poisson statistics were added to the data, the f(w^,w !f ) surface was determined 
from the noise-added data, and the reduced distribution was computed. This 
process was repeated until the spectrum of noise fluctuations in F(w lf ) no 
longer changed. With each computation of F(w n ) f the distribution was tested 
for a measurably significant positive slope. A total of 25 independent 
reduced distributions were computed from the noise-added data, and the bump on \ 

tail in Figure 2 and Figure 3 was present in 23 and 17 cases, respectively. 

We conclude that these are statistically significant bump-on-tail 
distributions. Furthermore, these reduced distributions have been found to be 
unstable according to the Penrose criterion [e.g., as in Krall and Tnvelpiece 
(1973)]- 


Discussion 

The magnetic field configuration and spacecraft location were particularly 
well suited for observing the foreshock boundary on November 6, 1977, the day 
of the measurements shown here. The geometry of the event shown m Figure 1 
remained fairly steady for several hours on that day, between 1000-1600 UT, 
approximately. With the leading edge of the foreshock near the spacecraft, 
short time scale changes m magnetic field direction caused the boundary 
pattern to sweep back and forth past the spacecraft many times during this 
period. Approximately 165 reduced distributions with a potential bump-on-tail 
were selected and 60 had measurable positive slopes that were significant to a 
68% level of confidence or better based on the error analysis outlined in the 
previous section. The remaining cases were rejected as resolved bump-on-tail 
distributions because the phase space sampling of the measurement and counting 
rate conspired to affect the reliability of the F(w Tt ) calculation near the 
bump of the distribution . 

The geometry of the boundary sketched in Figure 1 predicts that the 
distance between different parallel velocity characteristics increases with 
impact distance along the connecting field line to the shock and, therefore, 
the boundary in these circumstances should become more resolvable. During the 
time of these observations, the impact distance increased from 1 Re to 10 Re, 
and the number of cases in which two or more consecutive three second 
measurements 9 seconds apart in a given boundary crossing did in fact 
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increase. There was also a corresponding increase in consecutive measurements 
of bump-on-tail reduced distributions as the impact distance increased and a 
significant fraction of the total (22 of 60 cases) occurred m groups of two 
or more. This strongly suggests that the peaked reduced distributions are a 
persistent feature of the foreshock boundary. 

Conclusion 

The persistence of Penrose unstable beams on the reduced distribution 
function for macroscopic times (A t = 3-15 sec) > 300 u 1 is now an 

pe 

established fact at the electron foreshock. The kinetic explanation for the 
persistence of such beams remains for future theoretical work, since the 
simple oscillating two stream (OTS) scenario discussed by Filbert and Kellogg 
[197^] as adapted from Papadopoulos et al . , [1974] does not explain the 
stabilization of the beam at such low energies. Recent, more general OTS 
schemes [Freund et al. f 1981] of beam stabilization in the presence of 
suprathermal background solar wind electrons may prove to be more relevant to 
the actual circumstances of the foreshock. If this can be shown, ion-acoustic 
noise should be found in the locales of persistent beams. These and/or other 
explicit considerations of the inhomogeneous, driven character of the 
foreshock regime may lead to a kinetic explanation of the morphology presented 
above . 

Acknowledgments . We thank Richard Thompson for a computer routine which we 
used in transforming measurements of the velocity distr lbutions onto a 
rectangular grid suitable for contouring. 
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Figure Captions 

Sketch of electron foreshock boundary showing the spatial 
dispersion with respect to parallel velocity which generates the 
bump-on- tail reduced distributions. 

Contour map of the surface f(w x> w tl ) (upper panel) obtained from 
a three-dimensional measurement of the electron velocity 
distribution at the electron foreshock boundary on November 6, 
1977 at 15:39241-44 UT. The points underlying the contours 
locate the measurements in velocity space. The reduced 
distr ibution , F(w„) , (center panel) shows a bump on the tail at 

Q 

w„: -7 x 10 cm/s, and the vertical dashed lines outline the 
elemental strip of parallel velocity space producing the 

O 

positive slope. The enhancement in F at w |f 2 +7 x 10 cm/s, on 
the other hand, does not have a positive slope that can be 
supported by the propagation of errors analysis discussed m the 
text. (Lower panel) Contour map of w^ f(w x ,w M ) (normalized) 
showing the relative contribution of perpendicular velocities to 
F(w n ) at each parallel velocity. 

Another example of the bump-on-tail reduced distr ibutions 
measured at the foreshock boundary on November 6, 1977 at 
09:54:04-07 UT. The points underlying the contours locate the 
measurements in velocity space and show that there was good data 
coverage along the magnetic field direction. 
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